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Abstract. We propose that at the beginning of the Maunder minimum the poloidal field or 
amplitude of meridional circulation or both fell abruptly to low values. With this proposition, 
a flux transport dynamo model is able to reproduce various important aspects of the historical 
records of the Maunder minimum remarkably well. 
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1. Introduction 

One important aspect of the solar cycle is the Maunder minimum during 1645-1715 
when the solar activity was strongly reduced (Ribes & Nesme-Ribes 1993). It was not 
an artifact of few observations, but a real phenomenon (Hoyt & Schatten 1996). From 
the study of historical data (Ribes & Nesme-Ribes 1993), it has been confirmed that 
the sunspot numbers in both the hemisphere fell abruptly to nearly zero value at the 
beginning of the Maunder minimum, whereas a few sunspots appeared in the southern 
hemisphere during the last phase. It is also established from the cosmogenic isotopes 
data (Beer et al. 1998; Miyahara et al. 2004) that the cyclic oscillations of solar activity 
continued in the heliosphere at a weaker level during the Maunder minimum, but with a 
period of 13-15 years instead of the regular 11-year period. 

The most promising model of studying solar cycle at present is the flux transport 
dynamo model (Choudhuri et al. 1995; Durney 1995; Dikpati & Charbonneau 1999; 
Chatterjee et al. 2004). The main sources of irregularities in this model are the stochastic 
fluctuations in the Babcock-Leighton process of poloidal field generation (Choudhuri 
1992; Choudhuri et al. 2007) and the stochastic fluctuations of meridional circulation 
(hereafter MC) (Hathaway 1996). Therefore we propose that the polar field or amplitude 
of MC or both decreased at the beginning of Maunder minimum. With this proposition, 
we use a flux transport dynamo model to reproduce a Maunder minimum. The details 
of this work can be found in Choudhuri & Karak (2009) and Karak (2010). 

2. Methodology 

We cary out all the analyses with the flux transport dynamo model described in Chat- 
terjee et al. (2004). To reproduce the Maunder minimum, we perform the following three 
separate sets of experiments. Similar to Choudhuri et al. (2007), flrst, we decrease the 
polar fleld above O.8i?0 by a factor 7 after stopping the code at a solar minimum. We 
change the polar field by different amount in two hemispheres. In northern hemisphere, 
we take 7 = 0.0, whereas in southern hemisphere, it is 0.4. In addition, in this calculation. 
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Figure 1. Results covering the Maunder minimum episode, (a) The butterfly diagram, (b) The 
smoothed sunspot number. The dashed and dotted hues show the sunspot numbers in southern 
and northern hemispheres, whereas the sohd hue is the total sunspot number. From Choudhuri 
& Karak (2009). 

we decrease the toroidal field by multiplying it everywhere by 0.8 to stop the eruption for 
some time. This essentially reduces the strong overlap between two cycles in our model 
(see figure 13 of Chatterjee et al. 2004). After making these changes, we run the model 
for several cycles without any further change. In the second procedure for reproducing 
the Maunder minimum, we decrease the amplitude of MC vq abruptly to a very low 
value. After keeping it at low value for few years, we again increase it to the usual value 
but at different rates in two hemispheres. In the northern hemisphere, it is increased at 
slightly lower rate than the southern hemisphere. Note that in this case we have varied 
only Vq and no other parameters of the model. We have repeated this calculation in the 
low difFusivity model of Dikpati & Charbonneau (1999) too. Last, we have included the 
effect of the fluctuations of polar field along with the fluctuations of MC. We have run 
the model for different values of 7s from to 1 at each values of vq from a very low value 
to the average value. Then we find out the critical values of vq and the corresponding 7 
factor for which we get a Maunder-likc minimum. 



3. Results 

First, we discuss the results from the polar field reduction procedure. It is shown in 
Fig. [1] (see the caption also) . In order to facilitate the comparison with the observation 
data, we have marked the beginning of Fig. [I] to be the year 1640. From this figure we 
see that the sudden initiation but gradual recovery of Maunder minimum and the north- 
south asymmetry of sunspot numbers in the last last phase have been nicely reproduced. 



Maunder minimum 
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Figure 2. Parameter space of amplitude of MC (vo) and polar field reduction factor (7). The 
line shows the values of these parameters which are giving Maunder-like grand minima. 

We also find the cyclic oscillation of the poloidal field in the solar wind (shown in fig- 
ure 2 of Choudhuri & Karak 2009). This oscillation explains the cyclic behavior found in 
cosmogenic isotopes data. In this calculation, we have taken a = 21 m s~^ and turbulent 
difFusivity r/p of the poloidal field within the convection zone = 3.2 x 10^^ cm^ s~^. This 
combination of a and r]p gives the correct growth rate to produce Maunder minimum. 
In Choudhuri & Karak (2009), we list some other combinations which also reproduce 
maunder-like grand minima. 

In the second procedure of reducing vq alone, we again get results very similar to what 
are shown in Fig. [TJ if vq is abruptly made 10 m s~^ and then allowed to gain back its 
strength (see figures in Karak 2010). This procedure also reproduces all the important 
features of the Maunder minimum remarkably well. We underscore that in the advection- 
dominated model (e.g. Dikpati & Charbonneau 1999) we do not get this result. This is 
because the decrease of MC in the advection-dominated model produces more toroidal 
field and makes the cycle stronger, as analyzed by Yeates et al. (2008). 

We conclude that it is possible to reproduce the Maunder minimum by decreasing 
either polar field or vq to very low values. However, if we allow both the polar field and 
Vq to decrease simultaneously, then it seems possible to reproduce Maunder-like grand 
minima without making either the polar field or the MC so low. We have found that 
for each 7 there is a particular value of vq which can produce a Maunder-like minimum. 
The line in Fig. [2] shows the combinations of parameters giving a Maunder minimum of 
appropriate duration. Values of 7 and vq lying in the lower left of Fig.[2]give grand minima 
longer than the Maunder minimum, whereas values lying in the upper right give shorter 
minima. We should mention that Fig. [T]is reproduced just by taking 7 = 0.2 (actually 
7Ar = 0.0 and 75 — 0.4 is used) keeping vq unchanged. However, in this calculation, we 
have to reduce vq to around 13.5 m s~^ along with the polar field reduction to 0.2 in 
both hemispheres. This is because in earlier calculations we had reduced the toroidal 
field slightly along with the change of polar field. Here we have not done this reduction 
of the toroidal field. Additionally, here the values of a and rjp are slightly different giving 
different dynamo growth rate. 

4. Conclusion 

We have shown that most of the important features of the Maunder minimum can be 
reproduced quite well by assuming a simple ansatz that the polar field or the amplitude 



122 



Karak & Choudhuri 



of MC or both decreased significantly at the beginning of Maunder minimum. Because 
of our lack of knowledge about the physical conditions at the beginning of the Maunder 
minimum, we cannot say how exactly the Sun was driven to the Maunder minimum. 
However, we should mention that there are several independent studies (Wang & Sheeley 
2003; Miyahara et al. 2004; Passos & Lopeas 2009) suggesting that the amplitude of 
MC was weaker during the Maunder minimum. If this happens to be correct, then this 
study along with several other studies (Chatterjee et al. 2004; Chatterjee & Choudhuri 
2006; Jiang et al. 2007; Yeates et al. 2008; Goel & Choudhuri 2009; Karak & Choudhuri 
2011) indicate that the solar dynamo actually is diffusion-dominated and not advection- 
dominated. 
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